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The novel experimental realization of four-level optical quantum systems (ququarts) is presented. 
We exploit the polarization properties of the frequency non-degenerate biphoton field to obtain such 
systems. A simple method that does not rely on interferometer is used to generate and measure the 
sequence of states that can be used in quantum key distribution (QKD) protocol. 

PACS numbers: 42.50.-p, 42.50.Dv, 03.67.-a 



I. Introduction. Recently multi-dimensional (D>2) 
quantum systems or qudits attracted much attention in 
context of quantum information and communication. It 
is partly caused by fundamental aspects of quantum the- 
ory since the usage of qudits allows one to violate Bell- 
type inequalities longer than with two dimensional sys- 
tems (for the references see the review [if). Much inter- 
est in qudits also comes from the application point, espe- 
cially from the applied quantum key distribution (QKD). 
Multilevel systems are proved to be more robust against 
noise in the transmission channel, although measurement 
and preparation procedures of such systems seems to be 
much more technically complicated than in the case of 
qubits. Different aspects of the security of qudit-based 
protocols have been analyzed Q, Ll • Lately a proof-of- 
principle realization of a QKD protocol with entangled 
qutrits (D = 3) jj] and with qudits p| have been demon- 
strated. For the last years several elegant experiments 
have been performed where rHffer m]^ kinds of optical qu- 
dits were introduced 0, 0, S IM 0, EO, E3, Eal Recently 
an experiment that ensured the full co ntro l o v er a polar- 
ization qutrit state was demonstrated [12, [l3j . However 
polarization qutrits do not seem to serve as a practical 
candidate for the multilevel QKD, since it is impossible 
to achieve a demanded state using only SU2 transforma- 
tions that are done with linear optical elements. Mean- 
time transformations done by these elements are needed 
for practical realization of QKD protocols. 
In this Letter we present the results of the experimen- 
tal preparation, transformation, and measurement car- 
ried out with polarization based ququarts or quantum 
systems with dimensionality D = 4. 
II. Polarization ququarts and their properties. If 
one considers the two-photon field, generated via sponta- 
neous parametric down-conversion (SPDC) process, then 
the pure polarization state can be written as a superpo- 
sition of four basic components: 

I*} = c 1 \H 1 ,H 2 ) + c 2 \H l ,V 2 ) + c 3 \Vi,H 2 ) + c 4 \Vi,V 2 ). (1) 

Here Cj = jcjle 1 * 1 *,^ = 1,2,3,4) are complex probabil- 



verted photons. If the down converted photons have an 
only polarization degree of freedom, then a ququart state 
(JTJ converts to a qutrit state i.e. middle terms in be- 
come indistinguishable. In order to distinguish between 
these terms one must be able to distinguish between 
the down converted photons either in frequency, momen- 
tum or detection time. In experiments, described in this 
paper, we chose the collinear non-degenerate regime of 
SPDC, so twin photons that form a biphoton were hav- 
ing different frequencies and propagating simultaneously 
along the same direction. The sum of their frequencies 
was equal to the frequency of the pump, according to 
energy conservation. Polarization properties of this state 
can be described by Stokes parameters, which are defined 
as the mean values of Stokes operators, averaged over a 
state Q. Although the description of the light polariza- 
tion can be introduced only for the quasimonochromatic 
plane waves, it is possible to use P-quasispin formalism 
[14| to describe the polarization of arbitrary quantum 
beams with n modes, frequency or spatial. It is worthy 
to note that the frequency representation of the ququart 
l|TJl is isomorphic to the spatial one, when twin photons 
have the same frequencies but propagate in different di- 
rections ^3 . Frequently in the tasks of quantum commu- 
nication it is convenient to operate with states in a single 
spatial mode. For the two-frequency and single-spatial 
mode field, the formal definition of annihilation/creation 
operators is given by the sum of corresponding opera- 
tors in each mode. Also we take into account that these 
operators do commute for different frequency modes. So 
the Stokes parameters will contain time dependent terms 
exp{i(uj\— ui^jt) that describe "beats" of frequency modes 
and have no connection with the light polarization. How- 
ever, these terms vanish if one considers the finite detec- 
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\j, (J = 1,2) are the central wavelengths of down con- 
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tion time, that allows to classically average these " beats" : the matrix considered above has the form 



(So) = (a\a 



a 2 a 2 
b\bi 



b\h + b\b 2 ) - 2; 
(Si) = (a\ ai + a\a 2 - b\b± - b\b 2 ) = 2(| Cl | 2 - |c 4 | 2 ); 
(^2) = (a\bi + <4& 2 + b\a\ + bla 2 ) = 

2Re(c\(c 2 + c 3 ) + c 4 (c£ + eg)); 
(53) = (a\b\ + a\b 2 — b\a\ — b\a 2 ) = 
2Im(c\{c 2 + c 3 ) + c 4 (c* 2 + 4)). 

The polarization degree is given by 
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This definition of the polarization degree is just gener- 
alization of the commonly used classical one. It differs 
from the definition suggested in [lfj, where it serves as 
a witness of the state purity. In the case of polarization- 
based qutrit states 12, 13], the polarization degree P 3 
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2 - \c' 3 \ 2 + 2\c'*c' 2 + c'*c' 3 \ 2 with c[ = c^y/Zc'z = 
c 3 ,c 3 = C4 in Q was an invariant to unitary 
polarization transformations. Indeed it is impossible to 
prepare all demanded pure states, unless one uses inter- 
ferometric schemes with several nonlinear crystals |l2). 
In particular there is no way to transform the basic 
qutrit state |* 4 ) = \V, V) with P = 1 into the state 
l^i) = \H, V) with P — using retardation plates. How- 
ever, in the case of polarization ququarts, this quantity 
is no longer the invariant and can be changed by ap- 
plying local unitary transformations in each frequency 
mode. It can be achieved by using dichroic polarization 
transformers, which act separately on the photons with 
different frequencies. For example to transform the state 
|*4) = \Vx t ,V\ a ) into the state |* 2 ) = \H Xl ,V\ a ) one 
needs to use the retardation plate which serves as a half 
wave plate at Ai and as a wave plate at A2. The unitary 
transformation on the state is given by 4 x 4 matrix 
that is obtained by a direct product of two 2x2 matrices 
describing the transformation performed on each photon: 
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tu 
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with "transmission" and "reflection" coefficients 
ti = cos Si + i sin Si cos 2ai, = i sin Si sin 2ai, 
Si = 7r(n — n e )h/\i, where Si is its optical thickness, h is 
the geometrical thickness, and ai is the orientation angle 
between the optical axis of a wave plate and vertical di- 
rection. So for the unitary transformation of the ququart 
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The same dichroic plate performs direct and reverse 
unitary transformations between the polarization Bell 
states: |<i> + (~') and I'I'+M), which are represented by 
ququarts with 02=03=0, Ci=+(-)c4=;4= and Ci=c 4 =0, 
c 2 =+(-)c 3 =^ correspondingly. Similar transformations 
with frequency non-degenerate biphotons have been re- 
alized in . 

In general to prepare an arbitrary ququart state it is 
necessary to use four nonlinear crystals arranged in such 
a way that each crystal emits coherently one basic state in 
the same direction. But in particular cases a reduced set 
of crystals is quite sufficient to generate specific ququart 
states which can be used in applications. For example 
to prepare the polarization Bell states it was sufficient 
to use two crystals 0, . Moreover even single crystal 
allows one to prepare the useful subset of ququarts. 

III. Polarization ququarts in QKD protocol. 
The complete QKD protocol with four-dimensional po- 
larization states exploits five mutually unbiased bases 
with four states in each. In terms of biphoton states 
the first three bases consist of product polarization states 
of two photons while the last two bases consist of two- 
photon entangled states: 
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\VxV 2 ), 
\AiA 2 ), 
\LiL 2 ), 

\LiV 2 ); (6) 
\RiV 2 ), 
\ViL 2 ); 
\VxR 2 ). 



Here H,V, D, A, R, L indicate horizontal, vertical, +45 
and -45 linear, right- and left-circular polarization modes 
correspondingly, and lower indices numerate the fre- 
quency modes of two photons. It has been proved 
that using only first two or three bases is sufficient for 
the efficient QKD. Exploiting the incomplete set of bases 
one sacrifices the security but enhances the key genera- 
tion rate. That is why we will restrict ourselves to the 
first three bases and skip the rest two. As we will show 
experimentally, the set of twelve states can be prepared 
with a single non-linear crystal and local unitary trans- 
formations. We also present a measurement scheme that 
allows to discriminate the states belonging to one basis 
detcrministically thus allowing its implementation in re- 
alization of QKD protocol with polarization ququarts. 

IVa. Experimental setup. The experimental setup 
for generation and measurement of ququart states is 
shown at 
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FIG. 1: Setup for preparation and measurement of ququarts 

Fig-H The lOmW cw He-Cd laser operating at 325nm 
serves as a pump. A 15mm lithium-iodate type I crys- 
tal emits down-converted photons at central wavelengths 
Ai = 702nm, A 2 = 605nm within the spectral width 
of 2nm each, propagating collinearly with the pump, 
so the ququart state |Vi V2) is generated. Then, this 
state was subjected to transformations done by dichroic 
wave plate(s). Finally the state passed through the 
zero-order half- and quarter plates depending on which 
bases has been chosen. The measurement setup consists 
of a Brown-Twiss scheme with a non-polarizing beam- 
splitter. An interference filter centered at 702nm with 
a FWHM bandwidth of 3nm was placed in transmitted 
arm. This part of the scheme allows to reconstruct ar- 
bitrary polarization biphoton-ququart by registering co- 
incidence counts for different projections that are done 
by the polarization filters located in each arm 0] . Each 
filter consists of a zero-order quarter- and half-waveplate 
and a fixed analyzer. Two Si-APD's, linked to a coinci- 
dence scheme with 1.5 nsec time window, were used as 
single photon detectors. 

IVb. Experimental procedure. Let us consider as 
an example, the preparation of a state \H1V2) from the 
initial state IV1V2). This transformation can be achieved 
by a dichroic wave plate oriented at 45° that introduces a 
phase shift of 2tt between extra- and ordinary polarized 
photons at 605nm, and a phase shift of tt for the con- 
jugate photons at 702nm. Using quartz as birefringent 
material it is easy to calculate that the one of possible 
thicknesses of the wave plate that does this transforma- 
tion should be equal to 3.406mm j^. We used two plates 
PI and P2 with an effective thickness of 3.401mm. If 
then one can tilt these wave plates towards each other 
by a finite angle 9, then the optical thickness of the effec- 
tive wave plate will be changing and, at a certain value 
of 9, the desired transformation will be achieved. This 
corresponds to maximal coincidence rate when the mea- 
surement part is tuned to select \H1V2} state. Monitoring 
the coincidences, one can obtain the value of 9 for which 
the main maximum occurs. Then, fixing the tilting angle 
at this value, one can perform a complete quantum state 
tomography protocol in order to verify if the state really 
coincides with the ideal. In order to change the basis 



from I to II (III) , zero order half- (quarter) wave plates 
oriented at 22.5° (45°) were used. 

V. Results and discussion. Fig. [3 shows the coinci- 
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FIG. 2: Dependence of single counts (upper) and coincidences 
(lower) on tilting angle 6 

dences and single count rates versus the change of tilting 
angle 9 which determines optical thickness of the effective 
wave plate. If the measurement setup is tuned to select 
the state \H1~V2) then dependence of coincidences rate on 
the plate optical thicknesses Si is given by formula: 

Rcoin oc (ai t 6 2 t ai6 2 ) = sin 2 (<5i) cos 2 (c) 2 ), (7) 

whereas the single counts distribution in the upper chan- 
nel is given by 

hoinm oc (01*01) = sin 2 (<5i). (8) 

The solid lines at the Fig. show the theoretical curves. 
We performed tomography measurements for the main 
and additional maxima as well as for the minimum. The 
minima in coincidences occur when intensity in any chan- 
nel drops to zero, so it is not a necessary condition for 
distinguishing the orthogonal state to the one selected by 
given settings of polarization filters. Nevertheless accord- 
ing to calculations and our measurements the minimum 
in the coincidences at Fig. [21 exactly refers to the state 
\V\H 2 ). Starting from the | F702n.ro ^605nm) w e prepared 
and measured the whole set of the states from © belong- 
ing to the first three bases. The table I show components 
of the experimental (theoretical) density matrix as well 
as fidelity F defined by F = Tr(pthpexp) for some of the 
states. 

The main obstacle for the practical implementation 
of the free-space QKD protocol based on ququarts is 
that one needs to perform a fast polarization transfor- 
mation at the selected wavelengths. At the same time 
the method discussed in this Letter allows one to un- 
ambiguously distinguish all states forming chosen bases. 
The measurement set-up which has been already tested 
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State 




exp/ th\ 

P22 (P22) 




exp/ th\ 

P44 (Am) 


F 


\HxVz) 


0(0) 


0.984(1) 


0(0) 


0.016(0) 


0.98 


\V1H2) 


0.024(0) 


0.004(0) 


0.944(0) 


0.028(0) 


0.94 


\D 1 A 2 ) 


0(0) 


0.994(1) 


0(0) 


0.006(0) 


0.99 


|Ai£> 2 ) 


0.015(0) 


0(0) 


0.950(1) 


0.035(0) 


0.95 


\RiL 2 ) 


0(0) 


0.973(1) 


0.027(0) 


0(0) 


0.97 


|il^2> 


0.016(0) 


0(0) 


0.957(1) 


0.027(0) 


0.96 



TABLE I: Density matrix components and fidelities for se- 
lected transformations. 
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PBS 
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D4D2 


pu 


D4D1 


p22 


D3D2 


P33 


D 3 Di 


pa 





0.0(0) 


220 


0.973(1) 


6 


0.027(0) 





0.0(0) 



TABLE II: Coincidence rate and density matrix components 



of single photon detectors. 

To conclude we have suggested and tested a novel method 
of the preparation, and measurement of the subset of 
four-dimensional polarization quantum states. Since for 
this class of states the polarization degree is not invari- 
ant under SU2 transformations it is possible to switch 
between states using simple polarization elements. 
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FIG. 3: Measurement part at Bob's station 

in our experiments is shown on the Fig. |21 It consists 
of the dichroic mirror, separating the photons with the 
different wavelenghts, and a pair of polarization beam- 
splitters, separating photons with the orthogonal polar- 
izations. We would like to stress that using dichroic 
beam-splitter allows one to achieve 100% mode separa- 
tion efficiency which is not possible for qutrits. Four- 
input double-coincidence scheme linked to the outputs of 
single-photon detectors registers the biphotons-ququarts. 
For example for the first basis, the scheme works as fol- 
lows, provided that Bob's guess of the basis is correct: 
if state \H1H2) comes, then detectors D4, D2 will fire, 
if state \H1~V2) comes, then detectors D4, Dl will fire, 
if state \V1H2) comes, then detectors D3, D2 will fire, 
if state IV1V2) comes, then detectorsD3, Dl will fire. 
Same holds for any of the remaining correctly guessed 
bases, since the quarter- and half wave plates transform 
the polarization to HV basis in which the measurement is 
performed. Registered coincidence count from a certain 
pair of detectors contributes to corresponding diagonal 
component of the measured density matrix. So if the ba- 
sis is guessed correctly, then the registered coincidence 
count deterministically identifies the input state. We il- 
lustrate this statement by the table which shows total 
number of registered events per 30 sec for the input state 
\R1L2) measured in circular basis and calculated compo- 
nents of the experimental (theoretical) density matrix. 
Moreover, registering coincidences allows one to circum- 
vent the problem of the detection noise that is common 
for single-photon based protocols. If the coincidence win- 
dow is quite small, it is possible to assure a very low level 
of accidental coincidences for the usual dark count rate 
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